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1 .  INTRODUCTION 


As  part  of  the  DNA  stemming  and  containment  program  for 
underground  nuclear  testing,  Stanford  Research  Institute  (SRI) 
is  conducting  laboratory  investigations  to  develop  a  contain¬ 
ment  experiment  for  studying  residual  stress  fields  around 
exploded  cavities.  Systems,  Science  and  Software  (S1)  has 
been  asked  to  numerically  simulate  these  experiments  in  order 
to  increase  understanding  of  the  laboratory  results  and  to 
validate  our  capability  to  calculate  containment  related 
phenomena.  This  report  summarizes  one  part  of  our  calcula- 
tional  effort;  the  simulation  of  the  high  explosive  detonation 
and  the  subsequent  nonlinear  dynamic  processes  which  result  in 
the  formation  of  a  compressive  residual  stress  field  in  the 
grout  sphere  surrounding  the  exploded  cavity. 

The  experiment  is  shown  schematically  in  Figure  1  and 
the  experimental  apparatus  in  Figure  2  (Figures  1  and  2  have 
been  abstracted  from  an  SRI  report1) .  Step  A  of  the  experi¬ 
ment  involves  casting  a  12  inch  diameter  sphere  of  rock 
matching  grout  (RMG)  around  a  smaller  sphere  of  high  explo¬ 
sive  (PETN) .  The  explosive  and  detonator  are  sealed  off 
from  the  grout  by  a  thin  layer  of  lucite.  A  tube  containing 
the  detonator  wires  is  cast  into  the  grout  extending  almost 
to  the  anticipated  radius  of  the  exploded  PETN.  In  Step  B, 
the  grout  sphere  is  placed  in  the  tank  shown  in  Figure  2,  and 
a  constant  external  (overburden)  pressure  is  applied.  The 
water  in  the  tank  is  prevented  from  entering  the  grout  by  an 
impermeable  coating  on  its  surface.  In  Step  C,  the  high  ex¬ 
plosive  is  detonated  and  the  residual  cavity  gases  are  vented 
through  the  tube.  The  overburden  pressure  is  maintained 
during  the  drillback  to  the  cavity.  In  Step  D,  fluid  is 
pumped  into  the  cavity  at  a  constant  flow  rate  until  hydro¬ 
fracture.  The  major  laboratory  data  consists  of  graphs  of 
hydrofracture  pressure  vs.  volume  of  fluid  pumped  into  the 
cavity. 
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Figure  1.  Sequence  of  operations  in  containment  experiment. 


Hydrofracture  records  for  exploded  cavity  tests  art 
compared  with  records  for  unexploded  cavities,  i.e.,  for  qrout 
spheres  with  a  cast-in  cavity  and  access  tube.  The  earliest 
tests  showed  a  higher  pressure  was  needed  to  initiate  fracture 
for  an  exploded  sphere  than  for  an  unexploded  sphere  and  tha* 
the  shapes  of  the  records  after  fracture  initiation  were  quite 
different.  For  the  unexploded  case,  the  cracks  appeared  to 
propagate  at  lower  pressure  after  fracture  initiation  while 
for  the  exploded  spheres,  higher  pressures  were  recorded  after 
fracture  initiation.  These  test  results  appeared  to  confirm 
the  existence  of  a  compressive  residual  stress  field  around 
the  exploded  cavity.  However,  one  disturbing  note  was  tha+- 
while  the  exploded  sphere  usually  fractured  in  a  plane  not 
including  the  tube,  the  unexploded  sphere  usually  broke  in  a 
plane  including  this  tube.  In  addition,  the  test  results  were 
not  very  reproducible. 

In  the  last  two  years  major  improvements  and  modifica¬ 
tions  have  been  made  in  the  test  designs  including  varying  the 
amount  of  explosive  from  1/2  gm  to  .1/4  gm  to  the  present  3/8  gm, 
improving  the  charge  detonation,  reducing  the  lucite  shell  to 
the  present  11.4  mil  thickness,  changing  the  grout  mixture  to 
the  present  SRI  RMG  2C4 ,  improving  the  hydrofracture  fluxd 
measuring  system  to  reduce  air  bubbles,  and  eliminating  small 
inhomogeneities  in  the  unexploded  cavity  due  to  a  wire  mesh 
used  in  its  casting.  However,  the  basic  problem  of  reproduci¬ 
bility  of  the  exploded  sphere  results  remains.  In  addition 
the  latest  tests  show  very  little  difference  in  the  fracture 
initiation  pressures  for  exploded  and  unexploded  spheres,  which 
appears  to  cast  some  doubt  on  the  confirmation  of  the  residual 
stress  field,  the  primary  objective  of  the  experimental  study. 
However,  recent  tests  in  which  the  cavity  pressure  was  not 
vented  before  hydrofracture  result  in  higher  fracture  initiation 
pressures.  This  strongly  indicates  the  validity  of  the  resid¬ 
ual  stress  field  concept.  SRI  has  embarked  on  a  program  of 


6 


sensor  development  in  order  to  provide  additional  data  which 
may  be  compared  to  analytical  and  numerical  models  of  the  ex¬ 
periment.  These  sensors  include  a  strain  gauge,  a  residual 
stress  gauge,  a  surface  fracture  gauge,  a  stress  gauge  j.n  the 
wall  of  the  tank,  and  a  means  of  measuring  cavity  pressure 
before  venting. 

Here  we  present  the  results  of  our  one-dimensionai 
spherically  symmetric  dynamic  calculations  of  the  detonation 
of  the  PETN  explosive  and  the  formation  of  the  compressive 
residual  stress  field  in  the  grout  around  the  exploded  cavity. 
Results  are  presented  only  for  the  latest  test  configuration, 
a  3/8  gm  charge  of  PETN  contained  in  a  thin  lucite  shell,  sur¬ 
rounded  by  a  sphere  of  2C4  grout.  In  Section  2,  we  describe 
the  constitutive  models  and  material  properties  used  to  de¬ 
scribe  the  PETN  explosive,  the  lucite  shell,  the  2C4  grout, 
and  the  surrounding  water.  In  Section  3,  we  discuss  our  cal¬ 
culations  which  simulate  a  series  of  charge  calibration  experi¬ 
ments3  in  which  lucite  spheres  containing  3/8  gm  charges  of 
PETN  were  detonated  in  the  tank  filled  with  water.  Data  for 
these  tests  were  obtained  from  ytterbium  pressure  gauges  in 
the  water  surrounding  the  lucite  spheres  and  from  a  quartz 
gauge  mounted  in  the  tank  wall  directly  below  the  explosive. 
This  data  is  compared  with  our  calculations  and  with  other 
data  for  detonations  in  water. 

In  Section  4,  we  present  our  calculations  of  the  resid¬ 
ual  stress  field  in  the  grout  and  compare  our  results  with 
quartz  gauge  measurements  of  pressure  and  impulse  in  the  wall 
of  the  tank.  In  Section  5  we  discuss  a  simple  cal culational 
model  of  the  venting  and  hydrofracture  of  the  exploded  cavity 
which  allows  for  plastic  yielding.  This  model  assumes  a 
spherically  symmetric  hydrofracture  and  is  an  extension  of  the 
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analytical  work  of  K.  Narasimhan  for  linear  elastic  grout 
behavior.  (Two-dimensional  axisymmetric  hydrofracture  cal¬ 
culations  for  the  unexploded  sphere  case  have  been  discussed 


2.  CONSTITUTIVE  MODELS 


In  this  section  we  discuss  the  constitutive  models  and 
material  properties  used  to  simulate  the  behavior  of  the  PETN 
explosive,  the  lucite,  and  2C4  grout,  and  the  water  in  the  tank. 


2 . 1  PETN 

The  JWL  equation  of  state  developed  at  Lawrence  Livermore 
Laboratory^  was  used  to  describe  the  behavior  of  the  PETN  ex¬ 
plosive.  The  pressure  P  (Mb)  is  given  as  a  function  of  inter¬ 
nal  energy  E  (Mb  cc/cc)  and  V,  the  ratio  of  volume  to  initial 
volume  by 


v) 


-Rlv 


-r2v 

fc  + 


wE 

V 


where  A.  B,  ,  R 2  and  w  are  coefficients  of  the  fit. 

For  the  given  PETN  density  (pQ  =  1.0  gm/cc) ,  Finger7 
has  estimated  both  the  Chapman-Jouget  state  and  the  coeffi¬ 
cients  for  the  JWL  equation  of  state.  The  JWL  coefficients 
are 


A  =  2.372  Mb 
B  =  0.1061  Mb 
Rx  =  5.6 
R2  =  1.8 
w  =0.24 


The  Chapman-Jouget  state  is  given  by 

P  .  =  0.085  Mb 
C3 

D  =  0.555  cm/ys 

Eq  =  0.057  Mb  cc/cc 

T  =  2.69 

V  .  =  0.724 

c3 

where  D  is  the  detonation  velocity  and  (Ftl)/r  defines  the 
compression. 
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2.2 


LUCITE 


The  equation  of  state  used  to  describe  the  lucite  shell 

g 

was  developed  by  Duff  for  the  KO  code  at  Lawrence  Livermore 
Laboratory.  It  fits  shock  velocity-particle  velocity  data 
to  a  straight  line  of  the  form 

U  =  C  +  S  U 
s  p 

where  U  is  the  shock  velocity,  U  the  particle  velocity,  and 
s  p 

C  and  S  coefficients  of  the  fit  for  lucite  given  as 

C  =  0.251  cm/ us 
S  =  1.545 

The  pressure  is  calculated  from 

P  =  pQC  y  [l  +  u  +  S(S-l)u 

where 


and 

PQ  =  1.18  gm/cc 

Since  the  lucite  shell  was  sufficiently  close  to  the 
PETN  to  be  melted,  it  was  not  necessary  to  specify  its 
strength. 

2 . 3  GROUT 

Material  properties  tests  for  SRI  rock  matching  grout 
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(RMG  2C4)  were  conducted  by  Terra  Tek.  Physical  property 
and  ultrasonic  wave  velocity  measurements  were  made  as  well 
as  mechanical  property  tests  which  included  triaxial  compres¬ 
sion  and  uniaxial  strain  tests  to  0.5  Kbars  confining  pressure. 
Additionally  SRI  has  measured  unconfined  crush  strength.^- 
These  measurements  have  been  incorporated  into  our  equation 
of  state  for  grout. 
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Shear  Failure 


The  material  strength  Y  (the  maximum  stress  difference) 
is  chosen  to  be  a  function  of  P  as  described  in  Cherry,  et  al.^ 
and  is  given  by 

Y=Y  +  Y  /  2  -  )  P  <  P 

o  m  P  1  P  /  m 

m  \  m  / 

Y  =  Y  +  Y  P  >  P 

o  m  -  m 

where  Y  ,  Y  and  P  are  constants  of  the  fit  and  P  is  a 
o  m  m 

function  of  the  pressure  P  and  the  third  deviatoric  invariant 
(the  product  of  the  three  principal  deviatoric  stresses) . 

P  may  be  written  as 


In  addition,  the  stress  difference  Y  is  assumed  to  be 
zero  at  energies  above  the  melt  energy  (2xl010  ergs/gm)  and 
to  be  reduced  by  the  factor  1  -  e/em  for  energies  below  em- 

Based  on  both  the  uniaxial  strain  and  triaxial  test 

data,  the  maximum  value  of  Y  (YQ+Ym)  was  chosen  as  0.33  Kbars, 

and  P  ,  the  value  of  P  at  which  the  parabolic  strength  function 

is  replaced  by  Y  +Y  ,  was  taken  to  be  0.44  Kbars.  The  choice 
-1  o  m 

of  individual  values  for  Y_  and  Y  was  based  on  the  measured 

o  m 

unconfined  crush  strength  of  approximately  0.275  Kbars.  This 

is  equivalent  to  Y  of  0.21  Kbars  and  Y  of  0.12  Kbars  for 
^  o  m 

the  parabolic  form. 

Hooke's  law  is  used  to  obtain  an  initial  estimate  of 
the  stress  deviators.  Shear  failure  occurs  if  the  material 
strength  Y  evaluated  at  P  is  exceeded,  i.e.,  /3JJ  >  Y,  where 
is  the  second  deviatoric  stress  invariant  (one  half  the 
sum  of  the  squares  of  the  deviatoric  stresses) .  If  shear 
failure  occurs,  then  adjustment  of  the  stress  deviators  is 
required.  We  assume  that  each  deviatoric  stress  component  is 


reduced  by  multiplying  it  by  the  factor  f  given  by 


wh<  *e  b  =  dY/dp  is  the  slope  of  the  parabolic  strength  function. 

Tension  Failure 

Tension  failure  is  assumed  to  occur  in  an  element  if  a 
principal  stress  becomes  tensile  only  when  the  element  has 
previously  failed  in  shear.  We  then  apply  the  tension  failure 
model  proposed  by  Maenchen  and  Sack11  and  introduce  an  in¬ 
elastic  strain  normal  to  the  crack.  This  inelastic  strain  is 
just  sufficient  to  zero  the  tensile  stress  (the  failed  element 
cannot  support  any  tensile  stresses)  and  is  calculated  from 
Hooke's  law.  The  inelastic  strain  increments  are  accumulated 
during  each  cycle  in  which  the  principal  stress  is  calculated 
to  be  tensile,  thus  giving  an  estimate  of  the  crack  width. 

Once  this  unadjusted  principal  stress  becomes  compressive,  the 
crack  width  begins  to  decrease.  Crack  closure  continues  until 
the  inelastic  strain  becomes  zero.  When  this  state  is  achieved, 
the  element  is  able  to  support  a  compressive  stress. 


Equation  of  State  and  Pore  Crushup 

12 

The  Tillotson  equation  of  state  together  with  the  P-a 
porous  rock  crushup  model1^  were  used  to  describe  the  pressure 
of  a  rock  element  as  a  function  of  its  energy  and  density. 

The  pressure  of  the  material  without  air-filled  voids  is  given 
by  the  Tillotson  equation  of  state  as 


P  = 
s 


a  + 


V 


+  1 


ep  +  Ay  +  By 
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for  compressed  states  (p  >  o  ) ,  and  for  cold  expanded  states 
(p  <  pQ  and  e  <  es) .  The  form 


P 

v 


aep  + 


bep 


+  Aye 


-e2(v/vo-i) 


+  1 


Leon 


-6  (v/v  -1) 


is  used  for  expanded  states  (p  <  o  )  when  e  >  e'.  Here 

o  s 

VQ  =  l/pQ,  n  =  P/PQ  and  y  =  n-1.  A  simple  approximation  to 
the  phase  transition  from  liquid  to  vapor  is  obtained  for 
energy  states  eg  <  e  <  e^  when  o  <  pQ  by  calculating  the  pres¬ 
sure  as 

,e  -  eg)  Pv  +  (eg  -  e)  Pg  ] 

eg  is  the  specific  internal  energy  of  the  material  at  the 
vaporization  temperature  and  e^  includes  the  additional  heat 
of  vaporization  required  to  change  the  material  from  the 
liquid  to  the  vapor  state.  The  parameters  used  in  the  fit 
are  given  in  Table  1. 

For  a  material  containing  air-filled  voids,  the  pres¬ 
sure  may  be  obtained  from  the  pressure  given  by  the  Tillotson 
equation  of  state  through  the  use  of  the  P-a  crushup  model. 

The  distension  ratio  a  is  defined  as  the  ratio  of  the  density 
of  the  material  without  voids  to  its  density  with  voids  in¬ 
cluded.  Thus,  as  the  grout  is  loaded,  a  decreases  from  its 
initial  value  (l-<t>o)  where  <J)Q  is  the  air-filled  porosity, 
down  to  1.0  at  Pc,  the  pressure  at  which  all  air-filled 
porosity  is  irreversibly  removed.  Over  the  early  portion  of 
the  loading  (the  linearly  elastic  portion)  which  extends  up 
to  PQ,  the  elastic  pressure,  the  porosity  is  assumed  to  be 
completely  recoverable  upon  unloading.  Between  P  and  P_, 

6  C 

the  plastic  portion  of  the  crush  curve,  the  air-filled  porosity 
is  assumed  to  be  partially  recoverable  on  unloading. 


13 


Table  1.  Summary  of  2C4  Grout  Material  Properties  Data. 


Quantity 

Symbol 

Value 

Units 

Longitudinal  sound 
speed 

C 

3.27 

km/sec 

Shear  wave  speed 

Cs 

1.82 

km/ sec 

Density 

po 

2.16 

gms/cc 

Zero  pressure  bulk 
modulus 

K 

o 

135.57 

kbars 

Shear  modulus 

G 

71.55 

kbars 

Poisson's  ratio 

a 

0.275 

'  A 

150 

kbars 

B 

335 

kbars 

Bl 

5 

Coefficients  for 

$2 

5 

Tillotson  Equation  { 

.  a 

0.5 

of  State 

b 

0.633 

1  eo 

6 . OxlO11 

ergs/gm 

eSl 

3.5x1010 

ergs/gm 

1  X 

k  e 

s 

1.8X1011 

ergs/gm 

| 

f  Y 

1  xo 

0.21 

kbars 

Strength  Parameters  < 

1 

1  Pm 

0.12 

0.40 

kbars 

kbars 

1 

2xl010 

ergs/gm 

Air-filled  Porosity 

4>o 

0.007 

Elastic  Pressure 

Pe 

0.1 

kbars 

Crush  Pressure 

pc 

1.4 

kbars 
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Figure  3  shows  the  load-unload  curve  used  for  2C4  grout 
and  Table  1  gives  the  relevant  material  properties  data.  A 
table  of  values  for  a  as  a  function  of  specific  volume  V  may 
be  obtained  from  these  load-unload  curves  assuming  that  the 
unload  curve  represents  the  material  without  voids.  The 
curves  were  based  on  a  uniaxial  strain  test  performed  at 
Terra  Tek  for  2C4  grout,  in  which  the  sample  was  loaded  only 
up  to  a  pressure  of  0.5  kbars  and  then  unloaded.  It  appeared 
that  the  voids  were  not  completely  crushed  up  at  0.5  kbars. 
Fortunately,  the  same  sample  had  been  reloaded  up  to  4.0  kbars 
in  another  test.'*''3  Don  Gardiner  of  Terra  Tek  has  estimated 
the  total  air-filled  porosity  to  be  0.7  percent  by  adding  up 
the  voids  crushed  up  in  both  tests.  We  have  estimated  the 
crush  pressure  to  be  approximately  1.4  kbars  based  on  these 
tests.  The  coefficients  A  and  B  for  the  Tillotson  equation 
of  state  were  chosen  to  be  consistent  with  the  load-unload 
data. 

The  slope  of  the  elastic  portion  of  the  crush  curve 

was  chosen  so  that  the  zero  pressure  bulk  and  shear  moduli i 

are  consistent  with  the  measured  longitudinal  and  shear 

velocities  from  ultrasonic  tests.  The  choice  of  P  was  made 

e 

to  insure  elastic  behavior  at  the  overburden  pressure  of 
69  bars. 

2 . 4  WATER 

The  pressurized  water  in  the  tank  surrounding  the  grout 

sphere  was  modeled  using  an  analytical  equation  of  state  de- 

14 

veloped  by  Gurtman,  et  al.  which  treats  the  high  density 
region  near  the  Hugoniot  up  to  about  300  kbars.  This  poly¬ 
nomial  fit  cannot  however,  treat  expanded  states,  particularly 
if  phase  changes  occur.  These  restrictions  were  irrelevant  to 
the  grout  spheres  calculation  where  deviations  from  the  ap¬ 
plied  overburden  were  small  (less  than  one  overburden  pres¬ 
sure)  .  However,  for  the  charge  calibration  tests,  where  large 


.455  .460 

Specific  Volume  (cm3/?) 


Figure  3.  Load-unload  curves  for  2C4  grout 
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shock  pressures  were  expected  and  where  expanded  states  are 
possible  behind  the  shock  wave,  a  more  complete  equation  of 
state  for  water  may  be  desirable. 

The  AQUA  constitutive  relations  for  water, developed 
3 

at  S  ,  divide  the  pressure-volume  plane  into  several  regions. 
Within  the  steam  dome,  a  numerical  fit  to  data  given  in  the 
ASME  steam  tables  is  used.  In  the  region  bounded  by  the 
Hugoniot  the  200  kbar  release  adiabat,  and  the  liquid  satura¬ 
tion  line,  an  analytical  fit  is  used.  The  Hugoniot  is  de¬ 
scribed  by  another  analytical  fit  while  the  remainder  of  the 
P-V  plane  is  defined  by  large  data  tables  and  a  numerical 
table-lookup  procedure.  These  tables  are  also  extrapolated 
to  the  left  of  the  Hugoniot. 

Calculations  were  made  for  the  charge  calibration 
tests  using  both  equations  of  state  of  water  described  in 
this  section.  No  significant  differences  in  the  calculational 
results  were  obtained  when  the  analytical  equation  of  state 
was  replaced  by  the  AQUA  package.  Thus,  the  simpler  analytical 
equation  of  state  of  water  was  used  for  all  further  calcula¬ 
tions. 


3.  CHARGE  CALIBRATION  EXPERIMENTS 

A  series  of  experiments  were  made  to  calibrate  the  per¬ 
formance  of  the  3/8  gram  PETN  charges.  These  experiments  con¬ 
sisted  of  detonating  the  charge  configuration  shown  in  Fig¬ 
ure  4  in  water  in  the  tank  shown  in  Figure  2.  Pressure  histo¬ 
ries  at  various  locations  were  measured  by  means  of  a  piezo- 
resistive  ytterbium  gauge  and  a  piezoelectric  quartz  pressure 
transducer.  The  ytterbium  gauge  was  allowed  to  follow  the 
motion  of  the  water  so  that  a  direct  measurement  of  shock  wave 
pressure  was  obtained.  However,  the  quartz  gauge  was  rigidly 
mounted  in  the  bottom  of  the  tank  and  measured  reflected  pres¬ 
sure.  The  integrated  signal  is  therefore  reflected  impulse. 

Table  2  shows  the  results  of  12  charge  calibration  ex¬ 
periments  for  the  nominally  3/8  gm  PETN  charge  (the  actual 
weight  of  PETN  is  given  for  each  shot) .  Peak  incident  pres¬ 
sure  and  impulse  are  given  for  the  quartz  gauge  as  one-half 
the  corresponding  values  for  the  reflected  wave.  Also  pre¬ 
sented  in  Table  2  are  peak  incident  pressures  measured  by  the 
ytterbium  gauge  at  various  locations.  Note  that  the  quanti¬ 
ties  measured  by  the  two  gauges  do  not  have  a  consistent 
pattern  from  shot  to  shot.  For  example,  consider  shots  20, 

21,  and  22.  Based  on  charge  weight,  one  would  expect  the 
largest  peak  pressures  and  impulses  to  occur  for  shot  20. 
However,  the  maximum  impulse  at  the  quartz  gauge  is  approxi¬ 
mately  16  percent  less  than  for  shot  22.  (The  ytterbium  gauge 
gave  close  agreement  between  peak  pressures  for  these  2  shots.) 
Similarly,  a  comparison  between  shots  20  and  21  shows  slightly 
larger  impulse  at  the  quartz  gauge  for  shot  21,  but  28  percent 
smaller  peak  pressure  at  the  ytterbium  gauge.  The  variability 
in  these  results  seems  to  be  at  least  partly  due  to  ringing  in 
the  quartz  gauge.  Fortunately  most  of  the  integrated  impulses 
are  in  better  agreement  (shots  12  through  17) .  A  characteris¬ 
tic  of  the  data  is  that  all  records  integrated  electronically 
give  almost  the  same  impulse. 
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Table  2.  Summary  of  Results  for  3/8-Gram  Charges 


r 

Weight 

Shot 

of 

No. 

PETN 

(g) 

QUARTZ  GAUGE 


I  10 

i  11 


0.3630 

0.3644 

0.3634 

0.3566 

0.3586 

0.3672 

0.3587 

0.3619 

0.3701 

0.3679 

0.3674 

0.3794 


Distance 

from 

center  of 
charge  to 
gauge 
(in.) 


6-5/8 


YTTERBIUM  GAUGE 


Peak 

incident 

pressure 

(psi) 


3500 

4750 

3125 

4125 

4000 

4125 

2875 

3375 

3250 

4500 

4875 

4625 


50,000 (E) 
50, 500(E) 
50, 500(E) 
50,500 (E) 
50,500 (E) 
51,000  (E) 
53,500 (M) 
55,000 (M) 
62,000 (M) 
59,500 (M) 


7/16 

7/16 

7/16 

7/16 

7/16 

7/16 

7/16 

7/16 

1/4 

1/4 

1/4 

1/8 


54,100 

56,600 

55,900 

55,800 

50,750 

130,500 

102,600 

137,800 

144,700 


(E)  =  Electronic  integration 
Mechanical  integration 


E 


Two  calculations  have  been  performed  using  SKIPPER  to 
simulate  these  experiments.  One  used  the  AQUA  equation  oi 
state  for  water  and  the  other  a  simpler  analytical  fit.  Use 
of  the  different  water  equations  of  state  produced  only  negli¬ 
gible  differences  in  the  calculated  results.  The  constitutive 
modeling  used  for  the  PETN  and  lucite  have  been  discussed  in 
Section  2  of  this  report.  The  explosive  charge  configuration 
has  been  simplified  to  conform  to  the  spherical  geometry  of 
the  calculation.  A  3/8  gm  spherical  charge  was  used  in  the 
calculation  and  the  bridge  wire  was  eliminated.  This  gave  a 
radius  of  PETN  slightly  larger  than  reported  by  SRI  (0.447  cm 
rather  than  0.438  cm).  Since  the  outer  radius  of  the  lucite 
(0.476  cm)  was  the  same  as  in  Figure  4,  we  have  slightly  more 
PETN  (see  Table  2)  and  less  lucite  for  our  calculation. 

A  comparison  of  the  calculated  and  measured  peak  pres¬ 
sures  vs.  range  are  shown  in  Figure  5.  The  dashed  curve  indi¬ 
cates  the  relationship  predicted  for  explosive  charges  in 

1 6 

water  by  Langefors  and  Kehlstrom.  This  relationship  was 
compiled  from  many  sources  and  includes  measurements  over 
ranges  of  pressures  from  10  GPa  (100  kbars)  down  to  10  MPa 
(100  bars) .  Our  calculation  is  in  excellent  agreement  with 
these  predictions.  At  the  three  in-close  stations  (range  less 
than  2  cm) ,  the  calculated  peak  pressures  agree  well  with  the 
ytterbium  gauge  measurements.  The  slightly  low  measurement 
at  the  first  station  is  probably  not  significant  when  compared 
with  the  scatter  at  the  second  station. 

The  reflected  peak  pressures  measured  at  the  quartz 
gauge  are  considerably  higher  than  those  determined  by  the 
calculation  (approximately  a  factor  of  two  on  the  average) . 

A  portion  of  our  calculation  was  redone  with  finer  zoning. 

This  only  resulted  in  a  10  percent  increase  in  peak  reflected 
pressure  and  did  not  alter  the  total  impulse.  Figure  6  com¬ 
pares  the  pressure  history  calculated  at  the  quartz  gauge 
with  the  gauge  record  for  shot  15  of  Table  2.  The  gauge 
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Figure  5.  Calculated  peak  pressure  vs.  range  compared  to 
ytterbium  and  quartz  experimental  data  points. 

The  dashed  line  gives  the  relation  expected  for  an 
explosive  charge  in  water. 
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record  indicates  severe  nonphysical  ringing.  If  the  peaks 
and  valleys  are  averaged  and  a  smooth  curve  drawn  through  the 
average,  the  true  peak  pressure  is  probably  approximately 
30  MPa,  which  is  in  excellent  agreement  with  the  calculations. 
The  SKIPPER  calculation  shows  an  arrival  time  approximately 
9  microseconds  earlier  than  indicated  by  the  data.  This  may 
possibly  be  due  to  a  small  error  in  estimating  the  EDU  delay. 

If  the  calculated  pulse  were  displaced  9  microseconds  in  time, 
excellent  agreement  would  be  obtaired  between  the  calculated 
results  and  the  smooth  curve  drawn  through  the  ringing  quartz 
data. 

A  comparison  of  the  calculated  and  measured  impulse  at 
the  quartz  gauge  for  shot  15  (the  integral  of  the  waveform  of 
Figure  5)  is  shown  in  Figure  7.  Except  for  the  difference  in 
arrivals,  the  two  curves  are  in  remarkable  agreement,  partic¬ 
ularly  when  one  considers  the  shape  of  the  measured  pressure 
history,  and  the  fact  that  the  quartz  gauge  is  in  the  bottom 
of  a  cylindrical  tank,  whereas  the  calculation  assumes  spheri¬ 
cal  geometry.  At  times  greater  than  about  150  microseconds 
the  possible  influence  of  reverberations  off  the  sides  of  the 
cylindrical  tank  on  the  measured  signal  should  cause  our  cal¬ 
culated  impulse  to  differ  from  the  data.  The  calculation  was 
therefore  terminated. 

The  results  shown  in  this  section  have  shown  that  we 
can  accurately  calculate  the  detonation  of  the  PETN  and  the 
wave  propagation  through  the  lucite  and  water  using  the  SKIPPER 
code.  It  also  seems  clear  that  SRI  has  demonstrated  reproduci¬ 
bility  for  its  3/8  charges  at  least  to  within  20  percent  for 
the  measured  impulse.  The  measured  impulse  at  the  quartz 
gauge  should  be  a  useful  measurement  for  the  grout  spheres 
experiments  to  be  conducted  in  the  future  in  order  to  cali¬ 
brate  whether  the  variability  in  hydrofracture  records  is  due 
to  variations  in  the  PETN  detonation  and  dynamic  propagation 
through  the  grout. 
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Calculated  impulse  compared  to  integrated  quartz  pressure  history. 


4.  COMPUTATION  OF  THE  RESIDUAL  STRESS  FIELD 
IN  THE  GROUT 

Calculations  were  performed  of  the  detonation  of  the 
3/8  gm  PETN  charge  placed  in  a  12  inch  diameter  sphere  of  2C4 
grout.  Constitutive  models  and  material  properties  used  in 
the  calculations  have  been  presented  in  Section  2,  and  the 
configuration  of  the  charge  (PETN  plus  lucite)  has  been  dis¬ 
cussed  in  Section  3.  The  sphere  was  surrounded  by  water  at 
an  overburden  pressure  of  6.8  MPa,  and  a  rigid  boundary 
simulating  the  quartz  gauge  in  the  bottom  of  the  tank  was 
placed  17.8  cm  from  the  charge  center. 

The  calculations  showed  extensive  tensile  cracking  after 
passage  of  the  shock  wave  both  in  the  first  0.857  cm  (0.5  cm 
thick  after  cavity  expansion)  of  grout  adjacent  to  the  lucite 
and  in  an  interior  region  extending  radially  from  approximately 
5.]  to  7.5  cm.  These  cracks  were  all  in  the  radial  direction 
(the  hoop  stress  became  tensile) .  Those  cracks  near  the  ex¬ 
plosively  formed  cavity  closed  during  cavity  rebound.  How¬ 
ever,  the  permeability  of  this  region  almost  certainly  has 
been  increased  due  to  the  tensile  cracking.  Examination  of 
the  exploded  spheres  after  hydrofracture  shows  a  heavily 
stained  rubble  region  extending  approximately  1/4  inch 
(0.6  cm)  which  indicates  complete  penetration  by  the  hydro¬ 
fracture  fluid  that  could  be  due  to  increased  permeability. 

The  tensile  cracks  in  the  interior  region  of  the  grout 
do  not  close  completely.  Figure  8  shows  the  calculated  re¬ 
sidual  stress  field  in  the  grout.  The  radial  stress  ar  is 
continuous.  However,  the  plot  of  hoop  stress  shows  a 
small  bump  in  the  previously  cracked  region  at  a  radius  of 
approximately  6  cm,  just  before  the  remaining  open  portion 
of  the  crack  (where  ofc  is  zero) .  Both  stress  components  are 
equal  to  the  applied  overburden  in  the  water  at  15.2  cm. 
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(bars) 


Due  to  the  presence  of  a  rigid  boundary  m  t.,e  calcula¬ 
tions  at  the  location  of  the  quartz  gauge,  almost  constant- 
oscillations  with  variations  of  200  bars  in  peak  residual 
stress  were  calculated  as  a  wave  was  reflected  off  the 
boundary  propagated  inward  and  was  reflected  again  at  the 
cavity.  As  the  oscillation  passed,  the  remaining  tensile 
cracks  opened  and  closed.  These  oscillations  tended  to  con¬ 
tinue  throughout  the  grid  to  late  times.  To  obtain  our  final 
residual  stress  field,  the  coefficient  of  linear  artificial 
viscosity  was  increased  a  factor  of  10,  thus  damping  out  this 
oscillation . 

The  cavity  radius  corresponding  to  the  final  residual 
stress  state  was  calculated  as  1.21  cm  and  the  cavity  gas 
pressure  as  251  bars.  It  should  be  pointed  out  that  this 
calculated  cavity  pressure  is  considerably  larger  than  the 
average  fracture  initiation  pressure  for  all  3/8  gm  charge 
tests  which  is  185  bars  (2687  psi) .  We  believe  that  a  creep 
or  stress  relaxation  process  takes  place  in  the  time  between 
the  formation  of  the  residual  stress  field  and  the  hydro¬ 
fracture  of  the  grout  spheres.  This  process  reduces  the 
magnitude  of  the  compressive  residual  stress  field  so  that 
the  hydrofracture  can  occur  at  lower  pressures.  There  is 
some  experimental  evidence  to  support  this  hypothesis.  Grout 
spheres  that  are  left  standing  many  hours  before  hydrofractur¬ 
ing  give  lower  fracture  initiation  pressures.  In  addition, 
cavities  which  have  not  been  vented  (venting  should  accentuate 
the  creep  process)  give  significantly  larger  fracture  initia¬ 
tion  pressures. 

The  measured  cavity  radii  for  the  3/8  gm  charge  tests 
are  approximately  0.95  cm  (3/4  inch)  which  is  27  percent  lower 
than  the  calculated  value.  However,  the  radius  measurements 
take  place  after  the  sphere  is  cracked  open,  i.e.,  after  creep, 
venting,  and  hydrofracture  have  occurred.  There  is  little 
reason  to  believe  that  these  processes  have  not  altered  the 
cavity  radius  from  its  value  immediately  after  charge  detonation 
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One  purpose  of  the  quartz  gauge  is  to  test  and  validate 
predictions  in  order  to  gain  some  understanding  of  the  pheno¬ 
menology.  Figure  9  shows  the  calculated  overpressure  vs.  time 
at  the  quartz  gauge  and  Figure  10  the  calculated  total  impulse, 
obtained  by  integrating  the  data  of  Figure  9.  Recent  SRI 
data1^  for  tests  138  and  139  showed  peak  incident  overpressures 
of  530  (3.654)  and  550  psi  (3.792  MPa).  Maximum  incident  im¬ 
pulses  were  6200  (42.7)  and  6500  psi*ysec  (44.8  MPa*ysec).  The 
reflected  values  for  comparison  with  our  calculations  are 
twice  the  incident  giving  an  average  measured  peak  overpressure 
of  74.6  bars  (7.46  MPa)  and  an  average  maximum  impulse  of 
873  bars-ysec  (87.3  MPa-ysec) .  Our  calculated  maximum  values 
are  45.5  bars  overpressure  and  735  bars-ysec  maximum  impulse. 

As  discussed  in  Section  3,  the  pressure  measurement 
(see  Figure  6)  shows  a  great  deal  of  nonphysical  ringing.  Ex¬ 
trapolating  from  the  signal  shown  in  Section  3  (we  do  not  have 
the  pressure  histories  from  shots  138  and  139)  ,  the  reported 
peak  overpressure  could  be  high  by  almost  a  factor  of  2.  Thus, 
a  comparison  between  measured  and  calculated  peak  overpressures 
is  probably  not  meaningful  in  the  absence  of  a  full  pressure 
history.  However,  the  integrated  signal  (the  impulse)  has  been 
shown  in  Section  3  to  give  good  agreement  with  calculations 
even  for  a  ringing  pressure  history.  Therefore,  a  comparison 
between  measured  and  calculated  maximum  impulse  is  meaningful. 
Our  calculation  gave  a  15  percent  lower  impulse  than  the  aver¬ 
age  of  shots  138  and  139,  well  within  the  20  percent  variabil¬ 
ity  of  the  charge  calibration  test  results  given  in  Table  2  of 
Section  3.  Since  the  constitutive  models  and  materials  pro¬ 
perties  data  are  not  an  accurate  model  of  the  grout  to  better 
than  10  percent,  the  agreement  between  the  measured  and  cal¬ 
culated  impulses  appears  excellent. 
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Figure  10.  Impulse  at  the  quartz  gauge  (7")  for  3/8  gm  grout  sphere  calculation 


5.  A  NUMERICAL  MODEL  OF  THE  VENTING  AND  SPHERICALLY 
SYMMETRIC  FRACTURE  INITIATION 


A  majority  of  the  SRI  exploded  sphere  hydrofracture 
tests  show  one  major  crack  extending  out  of  the  cavity  region 
(as  indicated  by  the  dye  in  the  hydrofracture  fluid) .  There¬ 
fore,  a  numerical  model  of  these  tests  should  be  two-dimensional 
(see  Klein"’)  .  However,  some  of  the  SRI  tests  (shots  68  and  77) 
are  examples)  show  a  network  of  cracks  extending  radially  from 
the  cavity.  These  tests  also  require  considerably  greater 
pressures  in  order  to  hydrofracture  the  spheres.  In  order  to 
explain  these  tests,  a  simplified  one-dimensional  spherically 
symmetric  model  of  quasi-static  crack  propagation  through  a 

4 

residual  stress  field  was  developed  by  Narasimhan.  The  re¬ 
sults  obtained  from  the  model  provide  an  upper  bound  for  the 
hydrofracture  experiments  and  are  useful  in  describing  some 
of  the  qualitative  features  of  the  experiments,  i.e.,  the 
effects  of  impermeable  membranes  or  jackets  on  the  inside  or 
outside  of  the  sphere. 

4 

Several  of  the  assumptions  in  the  simplified  model 
are  in  conflict  with  the  residual  stress  fields  calculated 
here  using  the  SKIPPER  code  (see  Figure  8) .  In  the  simplified 
model,  the  residual  stress  field  for  most  cases  studied  is 
represented  by  an  analytical  expression  which  results  in  a 
zero  tangential  stress  at  the  outer  radius  of  the  sphere 
and  therefore  immediate  unstable  crack  propagation  for  the 
unjacketed  case.  A  few  cases  assumed  a  hoop  stress  equal  to 
one-half  the  overburden  pressure.  The  SKIPPER  calculations 
show  a  significant  hoop  stress  near  the  outside  of  the 
sphere.  Also,  the  simplified  model  assumed  the  grout  sphere 
to  be  uncracked  initially  before  hydrofracture,  but  the 
calculations  show  unclosed  cracks.  Finally,  excursions  from 
the  residual  stress  state  were  assumed  to  be  linearly  elastic. 
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Here  we  attempt  to  use  the  SKIPPER  code  to  simulate  the 
venting  of  the  high  pressure  explosive  gases  from  the  cavity 
and  the  subsequent  spherically  symmetric  hydrofracture  of  the 
grout  sphere  from  this  cavity.  In  this  way,  we  can  use  the 
actual  calculated  residual  stress  fields  including  the  tensile 
cracks  and  permit  the  grout  to  yield  plastically  during  hydro¬ 
fracture.  In  this  section,  results  are  presented  which  indi¬ 
cate  that  this  numerical  procedure  can  be  used  successfully  to 
model  hydrofracture.  We  should  emphasize  that  these  results 
are  only  preliminary.  Further  numerical  tests  are  presently 
being  made  to  determine  the  sensitivity  of  the  solutions  to 
zoning  and  to  unloading  and  reloading  rates. 

Figure  11  shows  the  hoop  stress  in  the  first  grid  zone 
of  the  grout  sphere  adjacent  to  the  cavity  as  the  pressure  in 
the  cavity  is  unloaded  (simulating  venting)  and  reloaded  (sim¬ 
ulating  the  onset  of  hydrofracture) .  Unload  was  accomplished 
by  reducing  the  cavity  energy  linearly  in  increments  of  time, 
thus  reducing  the  cavity  pressure  through  the  JWL  equation  of 
state  for  the  PETN  detonation  products.  Two  unloads  were  made 
down  to  a  pressure  of  1.0  bar,  one  in  a  time  of  30  micro¬ 
seconds  and  the  second  in  70  microseconds.  Point  A  in  Fig¬ 
ure  11  gives  the  initial  stress  state  in  zone  1  of  the  SKIPPER 
grid  after  the  dynamic  calculation  has  been  completed  (the 
cavity  pressure  is  approximately  the  same  as  the  radial  stress 
in  zone  1) .  As  the  cavity  pressure  is  decreased,  ofc  tends  to 
increase  (hoop  stress  becomes  more  compressive) .  Point  B 
indicates  the  onset  of  yielding  in  zone  1  for  the  70  micro¬ 
second  unload.  This  calculation  was  continued  for  approxi¬ 
mately  60  microseconds  after  the  unload  pressure  of  1  bar  was 
reached.  At  this  time,  the  hoop  stress  had  relaxed  down  to 
point  C  (afc  =  222  bars) .  Although  the  faster  unload  calcula¬ 
tion  (30  microseconds)  showed  a  different  unload  path,  the 
final  stress  state  was  in  fairly  good  agreement  (afc  = 

230  bars) . 
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In  these  preliminary  calculations,  no  attempt  was  made 
to  unload  in  a  time  comparable  to  the  actual  venting  time 
(tenths  of  seconds  to  seconds) .  Thus,  the  unload  was  not 
quasistatic  (the  transit  time  of  a  wave  from  the  cavity  to 
the  outside  of  the  sphere  and  back  is  approximately  100  micro¬ 
seconds  and  the  time  to  the  residual  stress  peaks  and  back  ap¬ 
proximately  60  microseconds) .  Negligible  changes  in  cavity 
radius  and  peak  residual  stress  resulted  from  the  unload.  We 
should  emphasize  that  time  dependent  creep  phenomena  which 
might  be  expected  to  reduce  the  stress  fields  were  not  included 
in  the  model.  It  remains  to  be  seen  what  the  effect  of  a  slower 
unload  would  be.  This  can  be  accomplished  by  a  rezone  of  the 
grid  near  the  cavity  which  would  greatly  increase  the  time  step 
(the  first  zone  in  the  preliminary  calculations  was  approxi¬ 
mately  0.005  cm  thick). 

Starting  from  the  stress  state  given  by  point  C  of 
Figure  11,  the  hydrofracture  of  the  grout  sphere  was  simulated 
by  a  time  dependent  pressure  boundary  condition  on  the  surface 
of  the  vented  cavity.  Loading  was  at  the  rate  of  72  bars  of 
cavity  pressure/msec,  i.e.,  much  more  rapidly  than  in  the  ex¬ 
periment.  There  was  an  initial  oscillation  in  in  zone  1  up 
to  pressures  of  40  bars  but  in  general  tended  to  decrease 
at  a  rate  one-half  as  great  as  the  loading  rate.  At  point  D 
along  the  loading  path  at  became  less  than  P  in  zone  1  for  the 
first  time.  At  this  point,  cavity  pressure  was  138  bars 
(2030  psi) ,  slightly  lower  than  measured  fracture  initiation 
pressures.  The  actual  fracture  initiation  pressure  depends 
on  inhomogeneities  in  the  grout  at  the  cavity  boundary  which 
will  initiate  cracking  once  a  stress  condition  is  reached 
which  will  propagate  a  crack.  It  has  been  speculated  that 
if  the  hydrofracture  fluid  (at  cavity  pressure)  can  find  a 
path  into  the  grout,  the  stress  state  given  by  point  D  in 
Figure  11  can  represent  the  minimum  value  of  the  fracture 
initiation  pressure.  Our  calculational  model  attached  no 


special  significance  to  this  stress  state  for  the  present 
calculations.  However,  the  cracked  region  adjacent  to  the 
cavity  described  earlier  for  the  dynamic  calculations  does 
indicate  a  path  into  the  grout. 

The  cavity  loading  was  continued  to  test  whether  the 

hoop  stress  a  became  tensile  (the  cracking  criteria  used  by 
4  t 

Narasimhan  ) .  Figure  11  shows  that  continued  to  decrease 
until  stress  state  E  was  reached,  but  that  from  then  on  a 
increased  as  loading  continued.  Point  E  indicates  the  onset 
of  yielding  in  zone  1.  At  this  point,  the  stress  path  changes 
from  a  path  characterized  by  constant  mean  stress  in  zone  1 
and  increasing  deviatoric  stress  in  the  radial  direction  (the 
reason  for  the  slope  of  one-half  in  Figure  11)  to  a  path 
characterized  by  increasing  mean  stress  as  well.  As  loading 
continues,  more  and  more  zones  begin  to  yield  and  ofc  continues 
to  increase.  The  cracking  criteria  of  Narasimhan  for  elastic 
loading  did  not  appear  to  be  reachable,  so  the  calculation 
was  terminated. 

If  future  sensitivity  tests  show  that  the  preliminary 
results  given  here  are  valid,  a  different  cracking  criterion 
must  be  developed,  possibly  based  on  some  critical  difference 
between  or  and  or  on  a  strain  criterion.  Once  a  cell  is 
cracked,  cavity  pressure  will  be  assumed  in  that  cell,  and 
the  calculation  continued.  We  hope  to  be  able  to  develop  this 
calculational  tool  to  investigate  the  grout  spheres  problem 
further.  Some  parametric  modification  of  the  residual  stress 
fields  before  unloading  to  simulate  creep  seems  also  worthy 
of  further  study. 

We  have  also  looked  at  the  possibility  of  investigating 
the  quasistatic  features  of  the  grout  spheres  experiments  with 
a  one-dimensional  finite  element  code.  At  present,  it  appc  rs 
more  cost  effective  to  use  the  SKIPPER  code  which  includes  the 
desired  constitutive  models  rather  than  to  incorporate  the* 
models  into  a  finite  element  code. 
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6. 


SUMMARY  AND  CONCLUSIONS 


In  this  report,  we  have  discussed  the  present  status  of 
one  aspect  of  our  numerical  simulation  of  the  SRI  grout  spheres 
containment  experiments.  We  have  demonstrated  our  ability  to 
correctly  calculate  containment-related  phenomena,  in  parti¬ 
cular,  a  high  explosive  detonation  and  the  subsequent  propaga¬ 
tion  of  nonlinear  waves  through  media  such  as  grout,  lucite, 
and  water.  Also,  we  have  introduced  a  simple  model  which  may 
be  useful  iri  explaining  some  features  of  the  experiment.  This 
model  uses  the  SKIPPER  code  to  examine  both  the  effects  of 
venting  on  the  calculated  residual  stress  fields,  and,  hope¬ 
fully,  also  to  study  the  onset  and  propagation  of  radial 
fractures  from  the  cavity  induced  by  the  hydrofracture  fluid. 

In  this  section,  we  summarize  our  results  and  conclu¬ 
sions  relative  to  the  various  phases  of  the  grout  spheres 
experiments. 

6.1  CHARGE  CALIBRATION  TESTS 

The  SRI  charge  calibration  tests  involved  exploding  the 
3/8  gm  PETN  charge  surrounded  by  a  thin  lucite  shell  in  a  water 
tank  and  measuring  peak  stresses  using  a  nearby  ytterbium  gauge. 
Stress  histories  were  also  measured  using  a  quartz  gauge  im¬ 
bedded  in  the  wall  of  the  tank.  Calculated  peak  stresses  vs. 
range  are  in  excellent  agreement  both  with  the  literature*^ 
and  with  the  ytterbium  data.  The  calculation  lies  well  within 
the  scatter  of  the  ytterbium  data.  Calculated  peak  reflected 
pressure  was  a  factor  of  2  lower  than  the  reported  data 
probably  due  to  severe  ringing  in  the  quartz  gauge.  If  the 
ringing  is  averaged  out,  the  data  then  is  in  excellent  agree¬ 
ment  with  the  calculated  history  (except  for  a  9  microsecond 
difference  in  signal  arrival  which  is  possibly  due  to  an  in¬ 
accurate  estimate  of  EDU  delay) .  The  integrated  signal  (the 
impulse)  is  in  excellent  agreement  with  the  calculated  impulse 
in  spite  of  the  ringing  gauge. 
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The  results  for  a  series  of  12  charge  calibration  tests 
indicate  that  the  measured  impulse  is  reproducible  to  within 
20  percent.  Thus,  it  will  be  very  useful  to  implace  a  quartz 
gauge  for  every  future  grout  spheres  tests  in  order  to  check 
out  the  charge  detonation. 

6.2  RESIDUAL  STRESSES  IN  A  GROUT  SPHERE 

Unfortunately,  very  little  data  exists  with  which  to 
compare  our  calculated  residual  stress  fields.  Our  calculated 
impulse  agrees  with  the  quartz  gauge  data  to  15  percent,  with¬ 
in  the  20  percent  variability  of  the  charge  calibration  test 
results.  As  before,  our  calculated  peak  pressures  are  much 
lower  than  the  reported  data  possibly  due  to  ringing  in  the 
quartz  gauge.  We  are  also  able  to  calculate  a  0.5  cm  thick 
region  of  tensile  cracks  adjacent  to  the  explosively  generated 
cavity  which  close  after  cavity  rebound.  These  cracks  result 
in  an  increased  permeability  which  is  seen  as  a  heavily 
stained  region  around  the  cavity  after  the  experiment  when 
the  spheres  are  split  open. 

The  calculated  cavity  pressure  before  venting  is  ap¬ 
proximately  30  percent  greater  than  the  average  fracture 
initiation  pressure  from  the  latest  SRI  exploded  sphere  tests. 
We  believe  that  a  creep  or  stress  relaxation  process  takes 
place  in  the  time  between  the  formation  of  the  residual  stress 
fie  Id  and  the  hydrofracture  of  the  grout  spheres  that  reduces 
the  magnitude  of  these  stresses  so  that  hydrofracture  can 
occur  at  lower  pressures.  Experiments  show  that  grout  spheres 
left  in  the  tank  for  hours  after  venting  hydrofracture  at  lower 
pressures  than  those  hydrofractured  immediately.  Also,  tests 
in  which  the  cavity  gases  are  not  vented,  and  in  which  reload 
takes  place  relatively  quickly  require  higher  fracture 
initiation  pressure.  Finally,  preliminary  results  using 
the  simple  calculational  model  described  in  Section  5  indi¬ 
cate  that  the  venting  itself  does  not  significantly  reduce 
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the  residual  stress  field.  These  all  indicate  the  presence 
of  a  time-related  stress  relaxation  process.  For  this 
reason,  we  strongly  recommend  that  all  future  explosive 
sphere  tests  use  unvented  cavities  and  that  they  begin  as 
soon  as  possible  after  detonation.  We  also  suggest  that 
creep  measurements  of  some  sort  be  made  in  order  to  guide 
development  of  a  creep  model. 
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